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The cell structure of clusters in the inner crust of a cold /3-equilibrium neutron star is studied 
within a Thomas Fermi approach and compared with other approaches which include shell effects. 
Relativistic nuclear models are considered. We conclude that the symmetry energy slope L may 
have quite dramatic effects on the cell structure if it is very large or small. Rod-like and slab-like 
pasta clusters have been obtained in all models except one with a large slope L. 



I. INTRODUCTION 

The inner crust of a neutron star lies between the neu- 
tron drip density {pd ~ 3 x 1CT 4 fm~ 3 ), denned as the 
density where the neutrons start to drip out from the 
nuclei of the crust, and the crust-core transition den- 
sity (p«8x 10~ 2 /m -3 ). In this region matter is most 
probably formed by a lattice of heavy and neutron rich 
nuclei immersed in a sea of superfluid neutrons and ul- 
trarelativistic electrons. Complex structures (e.g. rods 
and slabs) arc expected to be formed in the bottom part 
of the inner crust, the so-called pasta phase region, where 
the transition to the homogeneous core matter occurs. 

The first microscopic calculation of the inner crust 
structure was performed by Negclc and Vauthcrin in 
1973 [H. In this work, which is still used as a benchmark 
in neutron star calculations, the inner crust was studied 
in the Wigner-Scitz approximation, which divides the lat- 
tice in independent spherical cells each with the nucleus 
in the center surrounded by the electron and neutron 
gases. Nuclear matter was described with the Hartree- 
Fock (HF) approximation based on the Density Matrix 
Expansion [2J . The parameters of that model were ad- 
justed to reproduce the experimental binding energies of 
atomic nuclei and theoretical calculations of infinite neu- 
tron matter. The inner crust matter calculated in that 
work has for each cell a magic or semi-magic proton num- 
ber (i.e. Z = 40 and Z = 50), indicating that in these 
calculations there are strong proton shell effects, as in 
isolated atomic nuclei. 

Later Douchin and Haensel proposed a new study 
where the neutron star structure was calculated using the 
same Equation of State (EOS) in the whole star, from the 
outer crust to the core [3j. In that work nuclear matter 
was treated in the Wigncr-Seitz (WS) approximation us- 
ing the non-relativistic Compressible Liquid Drop (CLD) 
model and the SLy4 effective interaction J4|. The shell 
effects were neglected both for protons and neutrons. In 
the transition between the inner crust and the core five 
different pasta shapes were considered: spherical nuclei 
(droplets), cylinders of nuclear matter (rods), and plane 
slabs of nuclear matter (slabs) in a neutron gas, as well 
as cylindrical holes (tubes) and spherical holes (bubbles) 



in nuclear matter filled with a neutron gas. Within that 
model these authors found that in the whole inner crust 
the shape that minimizes the energy is the droplet and 
that those cells were characterized by a weak change in 
the proton number: from Zw40 near the neutron drip 
to Z ss 50 in the region of the transition to the core. 
Moreover because of the absence of proton shell effects 
in all the calculated configurations the proton number 
was not magic or semi-magic. 

The effect of pairing correlations on the structure of 
WS cells was investigated for the first time in [5[ within 
the Hartree-Fock-Bogoliubov approach (HFB). In the 
most recent version of these calculations these authors 
solved the HFB equations using an equation of state 
(EOS) mixture of the phcnomcnological functional of 
Fayans et al [fj, employed in the nuclear cluster region, 
and a microscopical functional derived from Bruckner- 
Hartree-Fock calculations in infinite neutron matter. In 
this framework it was found that the cells have not a 
magic or semi-magic number of protons and that pairing 
can change significantly the structure of the cells com- 
pared to HF calculations. 

The pairing effect was also studied by Gogelein and 
Mutter in 2007 Q- In this work these authors used a 
self-consistent Skyrme-Hartree-Fock (SHF) approach in 
order to study the cell structure in the pasta phase region 
(0.01 fm -3 < p < 0.1 fm~ 3 ), where three different pasta 
shapes were considered: droplet, rod and slab. The pair- 
ing correlations were evaluated within the BCS approach, 
assuming a density-dependent contact interaction, while, 
for the nuclear interaction, the SLy4 parametrization [4| 
was used. Contrary to what was obtained in [3[ in this 
work the calculations showed that all the three cell shapes 
appear in the inner crust before the transition to the 
homogeneous core. Moreover the calculated core-crust 
transition density was coherent with the value that was 
found in [8( for the SLy4 interaction within a dynamical 
spinodal calculation. 

Recently Grill et al. proposed a new inner crust struc- 
ture calculation based on the HFB approach J9|. This 
study was performed in the regions which are supposed 
to be formed by a lattice of spherical clusters. Thus the 
inner crust matter was divided in spherical cells treated 



in the Wigner-Seitz approximation and the structure of 
these cells (i.e. neutron number (N), proton number (Z) 
and cell radius (Rw S )) were obtained from the energy 
minimization at beta equilibrium. For the HFB calcu- 
lations it was considered a SLy4 interaction [4j in the 
particle-hole channel, while in the particle-particle chan- 
nel three zero-range density-dependent pairing forces of 
various intensities were used. With this model it was pos- 
sible to find very reliable results in the low density regions 
of the inner crust. Indeed the calculated structure is co- 
herent with the results of outer crust calculations in the 
literature (e.g. [101 )■ Moreover, in those regions a very 
weak dependence on the pairing interaction was found. 
On the other hand in the high density regions of the inner 
crust the imposed discretization of the free neutron gas 
generates an underestimation of the energy in the smaller 
cells: it was thus used an empirical correction |ll| which, 
however, was characterized by large fluctuations, so the 
structure of the cells in those regions was not guaranteed 
by those calculations. 

In this paper we present a calculation of the inner 
crust structure using a relativistic mean field density 
dependent Thomas- Fermi approach (TF) [lj|. A com- 
plete self consistent calculation is performed, namely, no 
paramctrization of the density and surface energy as in [7J 
has been used. This approach, which neglects the shell 
effect for both neutrons and protons, is less accurate than 
the HFB in the low density regions, where the nucleus in- 
fluences the cell properties. However, in the high density 
regions, when the free neutron gas becomes more impor- 
tant than the nucleus, the accuracy of the TF approx- 
imation can be higher than that of the HFB approach: 
indeed with this model the discretization of the free neu- 
tron gas is not necessary. Moreover, with the HFB ap- 
proach it is actually possible to consider just the spherical 
symmetry and so the droplet is the only available shape 
of the cells. On the other hand with the TF approach 
it is possible to study also the cylindrical and the plane 
symmetries, making available also the rod and slab cell 
shapes: we have performed our study also in the pasta 
phase region and determined the cell structure until the 
crust-core transition. 

The neutron inner crust is particulary sensitive to the 
density dependence of the symmetry energy [13 |. and, 
therefore, comparing results obtained with different nu- 
clear interactions will show how the symmetry energy 
affects the cell structure. This comparison is possible 
within a TF but would be prohibitive within a HFB cal- 
culation due to the excessive CPU time required. 



II. FORMALISM 

We will apply in the present study the self-consistent 
Thomas-Fermi calculation presented in Qj], [lj], within 
relativistic nuclear models with constant couplings and 
non-linear terms |15| . and with density dependent cou- 
plings [16}. In the relativistic mean- field theory protons 



and neutrons interact with and through an isoscalar- 
scalar field a, an isoscalar- vector field w M , an isovector- 
vector field p M and an isovector-scalar field S. Within the 
first class of models, that we will designate b y N on Linear 
Walecka Models (NLWM), we consider NL3 [T3| with non 
linear a terms, NL3 wp including also non-linear up terms 
which allow the modulation of the density dependence of 
the symmetry energy pj|, FSU [Hi and IU-FSU J2CJ with 
non-linear er, u and cop terms. These two parametriza- 
tions were constrained by the collective response of nuclei 
to the isoscalar monopole giant resonance (ISGMR) and 
the isovector dipole giant resonance (IVGDR). Within 
the second class of models with density dependent cou- 
plings we consider DD-ME2 [2l| and DD-ME5 [H]: DD- 
ME2, as all the non- linear parametrizations considered, 
does not include the 5 meson, and was adjusted to ex- 
perimental data based on finite nuclei properties; DD- 
ME<5 contains the S meson and was fitted to microscopic 
ab-initio calculations in nuclear matter and finite nuclei 
properties. 

Stellar matter will be described by a mixture of pro- 
tons, neutrons and electrons in chemical equilibrium. 
Electrons are described as a relativistic fermion gas which 
interacts with protons through the electromagnetic field 
A**. 

All the equations that allow the performance of the 
Thomas-Fermi calculation are derived from the La- 
grangian density 

C = ^ Ci + Ce +jC <x +£ u +£ p +jC S +£y+£nl, (1) 

where the nucleoli Lagrangian reads 
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iD* = id^-TM^-^T-p^-e^— —A", (3) 
M* = M-T a a-T S T-8, (4) 

and the electron Lagrangian is given by 

C e = i> e [ 7m (t0" + eA") - m e ] V> e . (5) 

The meson and electromagnetic Lagrangian densities are 
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where f^„ = dpUv-dvUp, R M! , = d ll p v -d v p lx -Y p {p lx x 
p v ) and F^y = d^A u — d^A^. The four coupling param- 
eters r a , T u , T p and Tg of the mesons to the nucleons 
are density dependent in the relativistic density depen- 
dent models considered, namely, DD-ME2 [2l| and DD- 
MEJ [221. The non-linear term C n i is absent in these 
models. In all the other models, NL3 [l7J], NL3.jp [Hj], 
FSU [H and IU-FSU 0, the couplings are constants 
and at least some of the non-linear terms of C n i are in- 
cluded. In the above Lagrangian density r is the isospin 
operator. 

The results obtained within the above relativistic 
mean-field models will be compared with the correspond- 
ing results calculated with the non-relativistic effective 
Skyrme interaction SLy4 [4| within two formalisms: a 
HFB calculation [<| and a CLD calculation 0). We 
will discuss how sensitive is the structure of the non- 
homogeneous inner-crust of a neutron star to the prop- 
erties of the EOS and the formalism used. In Table U 
the saturation nuclear matter properties and the a me- 
son mass are shown. The a meson mass has a strong 
influence on the nuclear surface energy and is included 
in the table to help the discussion. 

In Fig. [T] the symmetry energy and its slope L are 
plotted for all the models. The models considered have 
very similar values for the symmetry energy at satura- 
tion, namely, between 31.3 and 32.6 MeV except NL3 
that has a quite high value, 37.3 MeV. However, there is 
a larger dispersion of the symmetry energy slope L with 
values between 45 and 60 MeV, together with 118 MeV 
for NL3. All models behave in a similar way except NL3 
whose slope is much larger above p > po/i. SLy4 has the 
smallest slope L only at saturation density. FSU has the 
second largest slope only above the density p ~ 0.7 po- 
The properties of the pasta will reflect these facts, with 
IU-FSU and NL3 behaving in a quite different way, while 
all the other models showing similar results. The slope L 
has a direct influence on the surface energy tension and 
surface thickness of the clusters. 

A smaller L corresponds generally to a larg er surface 
tension and smaller neutron skin thickness [18| as can be 
confirmed comparing the surface tensions of the above 
models. In Fig. [5] the surface tension obtained from 
the derivatives of the meson fields, as indicated in [23j |. 
are plotted. The main differences between the models 
are due to the properties of the EOS at subsaturation 
densities: a smaller L dictates a larger surface tension 
in asymmetric matter: IU-FSU has the smallest L and 
largest surface tension; it is this effect that explains the 
difference between NL3 and NL3 wp . The dependence of 
the surface tension in symmetric nuclear matter on the 
model properties was well discussed in [24| , a smaller in- 
compressibility will give a smaller surface tension: this 
explains the small DD-ME5 value; a smaller value of m^ 
and a larger saturation density will give rise to a larger 
surface tension as in IU-FSU. 



III. RESULTS 

In our study of the inner crust structure we have 
considered three different cell shapes: the droplet, 
the rod and the slab. In most part of the in- 
ner crust, that we will designate by standard inner 
crust (p < 4.71 x 10~ 2 fm~ 3 ), we have limited our study 
just to the droplet configuration, so, as in [H, @, the lat- 
tice structure is described as a set of spherical cells, with 
radius R ws , treated in the Wigner-Seitz approach. The 
volume of a spherical cell with radius R ws is 
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However in the higher density regions, that we will desig- 
nate by pasta phase regions (p > 4.71 x 10~ 2 fm~ 3 ), we 
have taken into account all the three shapes. Rod cells 
have cylindrical shape with the radius R ws and length set 
for simplicity to I = 30 fm, which is ~ 1.5 Rws in region 
1 (the final results arc independent of this parameter) 



V r (R ws ) = InR 
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Slab cells have the shape of a parallelepiped with width 
and length set to I = 30 fm (the final results are indepen- 
dent on this parameter, too) and depth equal to 2R WSl 
so that 
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We consider these three different cell shapes at a given 
density and calculate the cell structure through the mini- 
mization of the energy per baryon under the condition of 
/3— equilibrium. For a fixed number of protons and neu- 
trons at a given density it is possible to obtain univocally 
the cell radius (R ws ) 
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and by the /3— equilibrium condition for each proton num- 
ber it is possible to calculate the neutron number. Fi- 
nally the cell structure is univocally defined searching 
the cell shape and the proton number that minimize the 
energy per baryon. In this minimization process we have 
treated Z and N as an integer (contrary to the approach 
of Douchin and Haensel 0]). The main effects of im- 
posing this condition are: smaller Wigner-Seitz cells and 
larger energies are obtained. 

In particular for each fixed density we have calcu- 
lated the cell shape, the proton and neutron number, 
the cell radius (R ws ), which is defined from the cell vol- 
ume ©, © and ©, the energy per baryon (E/A), the 
neutron chemical potential ((J.n) and the proton fraction 
{x = Z/{Z + N)). 



A. Standard inner crust 

The properties of the standard inner crust are reported 
in the TablcslHland lHll whcrc the proton number, neutron 



number, Wigner-Seitz radius, energy per baryon, neutron 
chemical potential and total proton fraction defined for 
twelve different values of the density defined in the sec- 
ond line of these tables. The first ten were introduced 
by Negele and Vautherin in [lj and we have added two 
more at low densities. In the following, we identify these 
densities by a label from 1 (high density close to the crust- 
core transition) to 12 (low density below neutron drip) 
and will refer to twelve density regions. In Fig. [3] and 
in Fig. @] we plot some of the properties of the Wigner- 
Seitz cells as a function of density. In particular, in Fig. 
[3] we show the neutron density at the cell center and 
the cell border, the cluster proton fraction at the cluster 
center and the neutron skin thickness = R n — i? p , with 
Rf = 5/3 < rf >, and in Fig. [4] we plot the Wigner-Seitz 
proton number Z, neutron number N, radius R ws , total 
proton fraction x, droplet number A and droplet proton 
fraction Z/A. The droplet nucleon number A has been 
estimated from the radius R ws and the neutron density at 
R ws , namely, A = Z + Ml - V i (R ws )p N (R ws )]. When- 
ever available data from [l], [j, Q have been included for 
comparison. In order to better understand the behaviour 
of the NL3 interaction we have added some points at high 
density (regions from 1 to 3). 

The main differences occur above p = 0.02 — 0.05 fm~ 3 . 
Below this density, except for the IU-FSU and NL3, all 
TF results lie between an upper and lower bound de- 
fined by the HFB and HF calculations. Above p = 0.05 
fin - non-spherical geometries may arise with lower en- 
ergy within the TF calculation, and this explains par- 
tially the differences. The CLD calculation for the SLy4 
interaction of [3j also follows the main trend obtained 
with the relativistic models. 

We will first discuss the proton number as a function 
of density. Different properties of the models explain 
the existing differences in the results. NL3 and IU-FSU 
have, respectively, the smallest and largest proton num- 
ber. This is a clear effect of the slope L, the large value of 
NL3 and the small value of IU-FSU, and corresponding 
surface tensions: in IU-FSU neutrons do not drip out so 
easily, the central droplet density is larger, and accord- 
ing to the liquid droplet model [25j, the proton number 
increases with the surface energy; the opposite occurs 
with NL3. NL3 wp only differs from NL3 in the isospin 
channel: a larger L clearly gives a larger proton number. 
The interaction DD-MEJ has, next to NL3, the smallest 
proton number and Wigner-Seitz radius in a large part 
of the density range considered, due to a smaller sur- 
face energy at large proton fraction which favors smaller 
droplets. FSU has a larger effective nucleon mass at sat- 
uration than all the other models, except IU-FSU. In [24[ 
it was shown that a larger effective mass favors a smaller 
surface thickness and, therefore, a larger nucleon number 
inside the droplet is expected, including a larger proton 
number. 

The results reported in Table [TT| show that in the TF 
calculations the proton number does not change consid- 
erably close to the transition to the outer crust and lies 



between 39 and 46. The proton number found with the 
HFB calculations converges to values around Z = 36, 
which is the proton number that characterizes the nuclei 
at the drip density (e.g. [HI)- This shell effect can not be 
reproduced by the TF calculations. The proton number 
obtained using the DD-ME2, DD-ME<5, NL3^ and FSU 
show a similar behaviour to the one found in [3j, where 
shell effects are also neglected: below p ~ 10 -3 fm~ all 
these models predict a larger proton number than the 
HFB [3] or the HF [l| calculations. This results are also 
similar to the ones calculated in [261 ] within the liquid 
drop model or in [27( within a TF calculation with en- 
ergy density functionals constructed in order to repro- 
duce nuclear matter properties. In this last work it is 
seen that Z is sensitive to the properties of the EOS. On 
the other hand in [28| . a relativistic Brueckncr-Hartree- 
Fock calculation was performed and quite small proton 
numbers were obtained, very similar to NL3. Our results 
also agree with the conclusions of [l3[ where, within a 
Thomas Fermi calculation applied to a macroscopic nu- 
clear model, it was shown that the proton number is 
larger for smaller values of L. 

In regions 1 and 2 the proton number found with the 
TF calculations has a consistent drop (Fig. 2]). This is 
a behaviour also occurring in other works (e.g. [lj, [5[, 
[281 ] and Q) and it is related with both the decrease of 
the Wigner-Seitz radius and the increase of the volume 
occupied with nuclear matter. 

The neutron numbers obtained with all interactions 
are very similar below p ~ 5 x 10 -3 fm and they are 
comparable to the ones calculated within the HFB Q 
and HF [l| calculations where shell effects are included, 
(sec Table [H| and Fig. [J). Above p ~ 5 x 10~ 3 fm~ 3 
model properties explain the existing differences. 

NL3 and DD-ME<5 interaction give neutron numbers 
smaller than the other interactions (except for NL3 in 
region 2). The neutron numbers found with NL3 wp (in 
the regions 1 and 2) and IU-FSU (in the region 1) are 
higher than what has been found with the other interac- 
tions. In fact, IU-FSU predicts the rod shape in region 1. 
If the droplet shape would have been considered in this 
region for IU-FSU the cells would have been character- 
ized by Z = 102, N = 2066 and R ws = 22.2 fm, much 
bigger than what has been found with all the other inter- 
actions. As a whole the neutron number calculated in all 
the considered interactions is coherent with the results 
found in the literature. It increases consistently with the 
density and only at very high density docs it decreases, as 
found in Q, 0] and [l[. We conclude that, as expected, 
the neutron number is not strongly linked to the shell 
effect. 

The TF and HFB neutron numbers are very similar 
just in the middle density range (corresponding to re- 
gions 5, 6 and 7). For larger densities (regions 2, 3 and 
4) the neutron numbers found with the HFB calculations 
are lower. Only in the region 3 the neutron numbers 
are similar, but this is a consequence of a big difference 
on the proton number. Finally, as expected, an impor- 



tant difference with respect to Q occurs in the low den- 
sity regions. As also found with the proton number, the 
HFB results converge to the closed shell neutron number 
(JV = 82) which characterizes the neutron d rip density 
and the initial part of the outer crust (e.g. [l0|), while 
the TF calculations, which does not include shell effects, 
cannot reproduce this behaviour. Furthermore, looking 
at the neutron chemical potential (Tab. HIT) ), we conclude 
that the neutron drip occurs above region 10 in the HFB 
calculations, due to the shell effects, while the TF cal- 
culations predict this transition in the region 12, which 
is characterized by quite low density. However, with the 
NL3 parametrization we have found a drip density similar 
to that found with the HFB calculations. TF results are 
very similar to those found by Negele and Vautherin [l|: 
only at the very high density (regions 1 and 2) there is a 
difference and the neutron number found in [l| is much 
larger than that found with TF. 

The cell neutron number is closely linked to the ra- 
dius of the Wigncr-Scitz cells (Tab. [ft] and Fig. gj}. The 
cell radius is quite independent on the nuclear interaction 
and the qualitative behaviour obtained in the TF calcu- 
lation is equal to what has been found with HFB [9| and 
HF [l[. Some quantitative differences identified are: a) 
the cells calculated with NL3 and DD-ME5 interaction 
are slightly smaller than those found with all the other 
interactions, which is coherent with the smaller neutron 
number predicted by these interactions (only region 2 is 
out of this trend in NL3); b) as with the neutron number 
the cell radius obtained using the TF and HFB approach 
are very similar in the middle density regions (from 5 to 
7) but differ in the low and high density ranges, where 
the HFB radius is smaller. 

In region 1, the lowest energy cell shape predicted by 
IU-FSU is the rod and, as it will be described in the 
next section, a change in the cell shape is always corre- 
lated with a drop in the value of the cell radius (Fig. [5]). 
However, this does not occur in this case because IU-FSU 
predicts values higher than those found with all the other 
interactions. 

The proton fraction calculated with all interactions is 
very similar (Tabs. IIII1 IVII and Fig. gj. Only the values 
calculated with the IU-FSU and NL3 interactions are dif- 
ferent from the general trend at higher densities (Tab. Hill 
in the regions from 1 to 4 and Tab. IVI[) . The behaviour 
of the proton fraction is similar to those found in the lit- 
erature [l|, @, 112, LU1 : its value decreases with increasing 
the density until, at high densities close to the crust-core 
transition, it changes its slope in order to reproduce the 
correct behaviour of homogeneous matter (Fig. @|. IU- 
FSU has the largest symmetry energy at p ~ 0.5po which 
favors isospin symmetry, NL3 has the smallest symmetry 
energy which favors a smaller proton fraction. In general 
a larger symmetry energy favors matter with a smaller 
isospin asymmetry. This is true for homogeneous matter, 
but the discussion of non-homogenous matter requires 
care: properties of clusters arc closely dependent on the 
surface energy which defines how favorable is the creation 



of clusters. 

The proton fractions obtained in [l2| and [7$ for the 
RMF model are smaller than the ones obtained in the 
present model, except for the NL3 parametrization. This 
is mainly due to the characteristics of the RMF used: 
both of them have a large symmetry energy slope at sat- 
uration and therefore, a smaller symmetry energy at sub- 
saturation densities, which favors large isospin asymme- 
tries. Similar results have been obtained in [29|. The cal- 
culated proton fraction in [7|] with the SLy4 parametriza- 
tion is consistent both with the CLD model and the HFB 
calculation, which have used the same parametrization, 
and are just slightly larger than the ones obtained with 
all the RMF models discussed in the present work except 
IU-FSU. In p} it was also shown that the shell effects 
were not very important and its main effect was to de- 
termine a slightly larger proton fraction. 

Some properties of the clusters reflect clearly the den- 
sity dependence of the symmetry energy (see Fig. [3] and 
H}. The cluster nucleon number above p = 5 x 1(U 3 fm -3 
is correlated with L while below that density it is not very 
sensitive to the model. The nucleon number of NL3 clus- 
ters decreases above 10" 2 /m~ 3 and at the crust-core 
transition is less than 100. On the other hand IU-FSU 
and SLy4, the models with the smallest slopes i, predict 
clusters with more than 500 nucleons close to the crust- 
core transition. All the other models predict clusters that 
do not go beyond 200-250 nucleons. 

The cluster proton fraction decreases, as expected, 
with density and estabilizes around 0.25 for models with 
L = 51 — 60 MeV. For both models with the smallest 
L, IU-FSU and SLy4 it decreases continuously reaching 
proton fractions below 0.15 while NL3 just before the 
crust-core transition has a strong increase of Z/A. We 
note, however, that as expected from the symmetry en- 
ergy close to saturation, NL3, with the largest symmetry 
energy, has the largest proton fraction at the center of the 
cluster, while IU-FSU has the smallest. The droplet over- 
all proton fraction is then defined by the clusters neutron 
skin (Fig. [3]), which is much larger in NL3. 



B. Pasta phase regions 

The properties of the pasta phase regions of the inner 
crust are reported in the tables from IIVI to I VIII and in 
Fig. [5] for the relativistic models. 

In the pasta phase regions the matter is studied con- 
sidering the three different cell shapes: the droplet, the 
rod and the slab. The transition densities between these 
different shapes are reported in table IVHI for all interac- 
tions. In this table we have also included the results of 
Ref. Q • All the three shapes appear in the inner crust 
except for NL3, which only predicts droplets. These re- 

13 and [U 



suits agree with those of Refs. [12j and [23[, where, us- 
ing the TF approach and the NL3 parametrization, it 
was predicted that in /?— equilibrium condition only the 
droplet cell shape would appear. Droplets were also the 



only configuration obtained in Refs. [3[ and [29|]. How- 
ever, contrary to |3| where the CLD was applied, in [7| all 
the three shapes were obtained with SLy4 using different 
frameworks, both TF and HF. In [13[ it was shown that 
models with a large L would not predict the existence of 
pasta shapes in /3-equilibrium matter. Thus we expect 
that the RMF parametrization used in [29[ has a large 
slope L. The IU-FSU shape transitions occur at densities 
similar to those obtained in pj for SLy4 within the HF 
calculation. These two models predict similar values of 
L. 

The crust-core transition densities found in this 
work, and indicated in Table IVII1 agree with the 
ones found in the literature within a dynamical spin- 
odal calculation: SLy4 at p t = 0.080 fm~ 3 , NL3 at 



p t = 0.054 fni , DD-ME2 at p t = 0.072 fnT J , FSU at 
p t = 0.074 fm~ 3 @, NL3 wp at p t = 0.0855 fm~ 3 [H and 
IU-FSU at p t = 0.087 fm -3 0. As expected results 
from the pasta calculation are just slightly larger. The 
crust core transition density is expected to occur within 
the metastable region between the spinodal and the bin- 
odal surfaces. However, for very asymmetric matter as 
/3-equilibrium matter, this region is almost nonexistent, 
since the binodal and the spinodal surfaces are very close. 
These results confirm the conclusions of reference |23| . 
namely that the estimation of the crust core transition 
from the dynamical spinodal is a good one. 

The proton and neutron numbers as well as the 
Wigncr-Seitz radius in the pasta phase are reported in the 
tables IIVI and [V] and in Fig. \5\ Comparing the behavior 
of the different models we conclude that: a) the NL3 and 
IU-FSU interactions are quite different from all the others 
with a much smaller or larger Z, N and R ws . The very 
different surface energy obtained within these two models 
explains this difference; b) all the other models behave 
in a very similar way, and differences may be explained 
comparing the symmetry energy slopes within the mod- 
els built using the same formalism, namely NLWM or 
density dependent hadronic models. A smaller L favors 
larger Z, N, and R ws in NL3 wp with respect to FSU. A 
smaller effective mass and incompressibility in DD-ME5 
favors a smaller surface energy and therefore, smaller Z, 
N, and R ws than in DD-ME2; c) the behavior of the sur- 
face energy with the proton asymmetry explains small 
differences between models. Close to the crust-core tran- 
sition asymmetries are large and models with larger sur- 
face energies will suffer shape transitions at larger densi- 
ties for similar cell sizes. 

At the shape transitions both the protons and neutron 
numbers have a sudden change because these numbers 
have a strong dependence on the cell volume (0 and 
(|8"|) and thus on the parameter I, which has been arbi- 
trary chosen. However this parameter has just an influ- 
ence on the proton and neutron number values, it cannot 
change their behaviour or the calculated matter proper- 
ties, namely their dependence on the density, the proton 
fraction, the energy per baryon or the chemical poten- 
tials. The cell radius decreases until densities close to 



a shape transition where it stabilizes or a small increase 
occurs just before the transition. The shape transition 
is then characterized by a decrease of cell size. Similar 
behavior was described in Qj], [29[ . 

A comparison with the results already published in the 
literature is not always possible. Indeed in Q just the 
standard inner crust has been studied, while the results 
found in |l| at very high density are not realistic because 
the droplet cell shape was imposed. In [3, Qjl, [29[ con- 
sistent calculations have been performed, but a detailed 
description of the cell structure (N,Z and R ws ) was not 
reported. The cell structure has been reported in Q 
(Tab. 1 and 2 of Ref. [3[]) where it is possible to see that 
the proton and neutron numbers have a behaviour similar 
to what has been found with our TF calculations: after 
an initial increase until quite high densities, then they 
decrease and finally they grow again until the crust-core 
transition (p = 7.5959 x 10 -2 fm -3 ). Also the cell radius 
has a behaviour very similar to the one obtained in this 
work: it monotonically decreases until densities close to 
the crust-core transition (p = 7.0154 x 10~ 2 fm~ ), then 
it grows until the transition. 

One interesting conclusion is that, except for NL3, all 
the models studied predict slab like configurations in /3- 
cquilibrium matter. According to [31] low energy collec- 
tive modes with an important contribution to the specific 
heat could be excited in these 'lasagna'-like pastas. 

The densit y in which non-spherical shapes appear was 
discussed in [25j, and an estimation of a nuclear filling 
fraction of 1/8 was obtained. If no dripped neutrons oc- 
cur this fraction translates into a average cell density of 
n,/8, where n, is the central density of the droplet. How- 
ever, this density will be larger if the dripped-neutron 
density is non zero. In fact, we confirm that the models 
with the lowest rod shape onset have the smallest neutron 
density at the cell border. This density is closely corre- 
lated with the slope L, with a smaller L corresponding to 
a smaller neutron density. NL3 is an exception because 
the crust-core transition occurs below this limit. The on- 
set of the slab geometry between different models follows 
a behavior similar to the rod onset. 



IV. CONCLUSIONS 

We have studied the inner crust properties of neutron 
stars within a self-consistent Thomas Fermi approach de- 
veloped in [12|,[lll for relativistic nuclear models. Several 
relativistic nuclear models have been used both with non 
linear meson terms and constant couplings, and with den- 
sity dependent coupling constants. The results have been 
compared with calculations obtained within the HFB, 
and the HF [l|, 0, [91 formalisms, with the compressible 
liquid drop model [3J] and with a macroscopic nuclear 
model [H. 

It has been shown that the main properties of the 
Wigncr-Seitz cells obtained within the HFB and HF 
formalisms are reproduced, namely, the average proton 



number and the neutron number and the Wigner-Seitz 
cell radius. As expected, proton shell effects are missing. 
The properties of the models used are reflected on the 
cluster structure. It was seen that a small symmetry en- 
ergy slope L gave rise to larger cells, with a larger proton 
and neutron number, while the opposite occurs for a large 
L. Models with a similar symmetry energy (~ 31 — 32 
MeV) and slope L (~ 50 — 60) at saturation density were 
shown to behave in a similar way, both in the droplet 
phase and the pasta phase regions. On the other hand 
models like NL3, with a very large symmetry energy and 
slope L and IU-FSU, with a quite small L, have shown 
quite different behaviors. NL3 did not present any pasta 
phases in the inner crust of /^-equilibrium matter, and 
predicted the smallest proton and neutron numbers, and 
Wigner-Seitz radius in almost all the inner crust range 
of densities. On the other hand, IU-FSU predicts a quite 
low density for the onset of the pasta phase, where all 
the other models still predict the existence of droplets. 
The occurrence of the slab shape occurs at a lower den- 
sities than in all the other models. However, the IU-FSU 
crust-core transition density is the largest one and above 
0.01 fm -3 IU-FSU presents the largest clusters with more 
than the double of nuclcons. All the models, except 
NL3, predict the existence of slab like configurations in j3- 
cquilibrium matter. These 'lasagna'-like structures may 



have an important contribution to the specific heat of the 
crust [3l[ . 

The size and composition of the clusters will have an 
im por tant effect on the transport properties of the crust. 
In [32j it was shown how the pasta structures could affect 
the neutrino transport, namely a more uniform distribu- 
tion of matter, as occurs in NL3 with a larger neutron 
drip, or a larger range with non-spherical pasta struc- 
tures could reduce the cross section for elastic neutrino 
scattering from pasta phases via weak neutral current, 
and therefore the neutrino opacity. 

The effect of the pairing correlations, which were miss- 
ing in the present work, on the inner crust clusters and 
the effect of the size and composition of the clusters on 
the transport properties of the crust will be investigated. 
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Po [fm- 3 ] 


E [MeV] 


K [MeV] 


E sym [MeV] 


L [MeV] 


M*/M 


m a [MeV] 


SLy4 


0.159 


-15.97 


229.8 


31.8 


45.3 


0.695 


- 


NL3 


0.148 


-16.24 


270.7 


37.3 


118.3 


0.600 


508.2 


DD-ME2 


0.152 


-16.14 


250.8 


32.3 


51.4 


0.609 


550.0 


DD-ME5 


0.152 


-16.12 


219.1 


32.4 


52.9 


0.572 


566.2 


NL3„ P 


0.148 


-16.30 


272.0 


31.7 


55.2 


0.600 


508.2 


FSU 


0.148 


-16.30 


230.0 


32.6 


60.5 


0.620 


491.5 


IU-FSU 


0.155 


-16.40 


231.2 


31.3 


47.2 


0.620 


491.5 



TABLE I: Nuclear matter properties at the saturation density (density, binding energy, incompressibility, symmetry energy, 
symmetry energy slope and effective mass) and the a meson mass. 



Regions 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


p [10- 3 fm" 3 ] 


47.1 


20.2 


8.83 


5.72 


3.70 


1.58 


0.871 


0.595 


0.396 


0.276 


0.259 


0.188 


Z 
HFB 




40 


54 


40 


46 


50 


50 


36 


38 


36 


38 


38 


N&V 


40 


50 


50 


50 


50 


40 


40 


40 


40 


40 


- 


- 


NL3 


16 


27 


32 


34 


36 


38 


38 


39 


39 


39 


39 


39 


DD-ME2 


33 


39 


42 


43 


43 


43 


43 


43 


43 


42 


43 


42 


DD-ME5 


32 


39 


42 


42 


42 


42 


41 


41 


41 


41 


41 


40 


NL3„ P 


39 


42 


44 


44 


44 


44 


44 


44 


43 


43 


43 


42 


FSU 


37 


44 


46 


46 


46 


46 


46 


46 


45 


45 


45 


44 


IU-FSU 
N 
HFB 


82 r 


63 


54 


52 


50 


48 


47 


47 


47 


46 


46 


45 


— 


1018 


1324 


732 


740 


454 


316 


174 


120 


<S2 


82 


82 


N&V 


1460 


1750 


1300 


1050 


900 


460 


280 


210 


160 


140 


- 


- 


NL3 


1019 


1816 


1202 


925 


718 


410 


259 


193 


134 


98 


96 


89 


DD-ME2 


1116 


1463 


1213 


1016 


808 


482 


313 


229 


160 


110 


106 


92 


DD-ME5 


1107 


1301 


1087 


908 


739 


448 


286 


211 


147 


105 


99 


86 


NL3„ P 


1300 


1590 


1291 


1045 


821 


476 


308 


224 


152 


108 


102 


92 


FSU 


1192 


1482 


1243 


1035 


835 


503 


328 


241 


164 


116 


110 


96 


IU-FSU 


1655 r 


1417 


1188 


1031 


840 


510 


332 


247 


173 


120 


113 


95 


Rw* [fm] 


























HFB 


- 


23.2 


33.4 


31.8 


37.0 


42.4 


46.4 


43.8 


45.6 


46.8 


48.0 


53.4 


N&V 


19.6 


27.6 


33.1 


35.7 


39.3 


42.2 


44.3 


46.3 


49.2 


53.6 


- 


- 


NL3 


17.4 


27.9 


32.2 


34.2 


36.5 


40.8 


43.3 


45.3 


47.1 


49.1 


49.9 


54.6 


DD-ME2 


18.0 


26.1 


32.4 


35.4 


38.0 


43.0 


46.0 


47.8 


49.6 


50.8 


51.6 


55.5 


DD-ME5 


17.9 


25.1 


31.3 


34.1 


36.9 


42.0 


44.8 


46.6 


48.4 


50.1 


50.5 


54.3 


NL3„ P 


18.9 


26.8 


33.0 


35.7 


38.2 


42.9 


45.9 


47.6 


49.0 


50.7 


51.1 


55.5 


FSU 


18.4 


26.2 


32.7 


35.6 


38.5 


43.7 


46.8 


48.7 


50.1 


51.8 


52.3 


56.3 


IU-FSU 


19.8 r 


26.0 


32.3 


35.6 


38.6 


43.9 


47.0 


49.1 


51.0 


52.3 


52.7 


56.3 



TABLE II: Cell proton number, neutron number and Wigner-Seitz radius for different nuclear interactions at the considered 
densities. HFB and N&V refer respectively to the results found in [9| and [J]. In the region 1, for the IU-FSU interaction, the 
cell shape is the rod. 
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Regions 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


p [l(T 3 fm- 3 ] 


47.1 


20.2 


8.83 


5.72 


3.70 


1.58 


0.871 


0.595 


0.396 


0.276 


0.259 


0.188 


E/A [MeV] 


























HFB 


- 


4.747 3.025 


2.283 


1.645 


0.612 


-0.051 


-0.490 


-1.063 


-1.691 


-1.830 


-2.441 


N&V 


6.428 


4.097 2.610 


1.996 


1.465 


0.541 


-0.050 


-0.462 


-0.962 


-1.425 


- 


- 


NL3 


4.491 


3.054 


2.438 


2.096 


1.735 


0.993 


0.422 


-0.004 


-0.545 


-1.130 


-1.244 


-1.801 


DD-ME2 


6.847 


4.494 


3.130 


2.560 


2.051 


1.174 


0.584 


0.170 


-0.341 


-0.890 


-1.002 


-1.581 


DD-ME5 


6.986 


4.796 


3.313 


2.681 


2.123 


1.186 


0.569 


0.139 


-0.391 


-0.959 


-1.075 


-1.666 


NL3„ P 


6.750 


4.354 


3.020 


2.474 


1.984 


1.122 


0.526 


0.100 


-0.430 


-1.003 


-1.119 


-1.702 


FSU 


6.987 


4.660 


3.202 


2.601 


2.072 


1.173 


0.575 


0.155 


-0.362 


-0.917 


-1.030 


-1.609 


IU-FSU 


8.531 


5.630 


3.645 


2.884 


2.252 


1.250 


0.623 


0.194 


-0.325 


-0.881 


-0.994 


-1.587 


Pn [MeV] 


























HFB 


- 


7.178 


4.600 


3.713 


3.117 


1.730 


1.036 


0.724 


0.327 


-2.660 


-3.278 


-3.287 


N&V 


10.900 


6.500 


4.200 


3.300 


2.600 


1.400 


1.000 


0.600 


0.300 


0.200 


- 


- 


NL3 


7.285 


3.242 


2.550 


2.261 


1.932 


1.253 


0.816 


0.556 


0.288 


-0.031 


-0.189 


-0.809 


DD-ME2 


10.118 


5.842 


3.883 


3.147 2.511 


1.509 


0.971 


0.677 


0.398 


0.145 


0.094 


-0.442 


DD-ME5 


9.717 


6.218 


4.224 


3.394 


2.679 


1.572 


0.995 


0.689 


0.397 


0.144 


0.094 


-0.515 


NL3u,p 


10.050 


5.712 


3.715 


3.009 


2.409 


1.455 


0.935 


0.646 


0.365 


0.105 


0.046 


-0.568 


FSU 


9.938 


6.130 


4.040 


3.244 


2.569 


1.526 


0.974 


0.677 


0.390 


0.133 


0.086 


-0.526 


IU-FSU 


11.659 


7.907 


4.933 


3.822 


2.934 


1.671 


1.051 


0.731 


0.432 


0.179 


0.131 


-0.425 


HFB 


— 


0.038 


0.039 


0.052 


0.059 


0.099 


0.137 


0.171 


0.241 


0.305 


0.317 


0.317 


N&V 


0.027 


0.028 0.037 0.045 


0.053 


0.080 


0.125 


0.160 


0.200 


0.222 


- 


- 


NL3 


0.015 


0.015 


0.026 


0.035 


0.048 


0.085 


0.128 


0.168 


0.225 


0.285 


0.289 


0.305 


DD-ME2 


0.029 


0.026 


0.033 


0.041 


0.051 


0.082 


0.121 


0.158 


0.212 


0.276 


0.289 


0.313 


DD-ME5 


0.028 


0.029 


0.037 


0.044 


0.054 


0.086 


0.125 


0.163 


0.218 


0.281 


0.293 


0.317 


NL3 W/3 


0.029 


0.026 


0.033 


0.040 


0.051 


0.085 


0.125 


0.164 


0.221 


0.285 


0.297 


0.313 


FSU 


0.030 


0.029 


0.036 


0.043 


0.052 


0.084 


0.123 


0.160 


0.215 


0.280 


0.290 


0.314 


IU-FSU 


0.047 


0.043 


0.043 


0.048 


0.056 


0.086 


0.124 


0.160 


0.214 


0.277 


0.289 


0.321 



TABLE III: Cell energy per baryon, neutron chemical potential and proton fraction for several nuclear interactions at the 
considered densities. HFB and N&V refer respectively to the results found in [f| and [J|. In the region 1, for the IU-FSU 
interaction, the cell shape is the rod. 
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FIG. 1: (Colors online) Symmetry energy E sym (a) and its slope L (b) as a function of the density p. 
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p 


- 3 ] 






Z 












N 








[fm 


NL3 DD-ME2 DD-ME5 NL3„ P 


FSU IU-FSU 


NL3 


DD-ME2 DD-ME5 NL3^ P 


FSU IU-FSU 


8.94 x 


10" 2 
10 -2 








— 




- 








— 




- 


8.88 x 


173 


3594 


8.35 x 


10" 2 
10" 2 








- 




117 
114 








- 




2418 


8.29 x 


67 


1680 


2355 


7.82 x 


10~ 2 


- 


- 


- 


60 


- 


101 


- 


- 


- 


1557 


- 


2071 


7.66 x 


10" 2 
10 -2 
10~ 2 
10" 2 
10" 2 


- 


- 


- 


59 
59 


— 


99 
97 
95 
94 
93 


- 


- 


1396 
1385 
1373 
1095 


1548 
1565 
1319 
1280 
1267 


1503 
1439 
1401 


2027 


7.51 x 


49 

48 
47 


1981 


7.35 x 


49 
47 
46 


56 
53 

51 


1935 


7.20 x 


54 
51 


1467 
1403 


1910 


7.04 x 


37 


1885 


6.88 x 


10~ 2 
10~ 2 


— 


50 


36 
35 


45 

45 


50 

50 


92 
92 


— 


1394 


1079 
1061 


1254 
1269 


1391 

1407 


1861 


6.73 x 


41 


1158 


1856 


6.57 x 


10~ 2 


- 


39 


35 


44 


41 


91 


- 


1116 


1074 


1255 


1170 


1833 


6.42 x 


10~ 2 


- 


38 


35 


44 


40 


91 


- 


1102 


1086 


1271 


1153 


1830 


6.26 x 


10~ 2 
10~ 2 
10~ 2 


- 


38 
38 


35 


42 
41 

40 


40 
39 
39 


90 
99 

96 


- 


1115 
1131 


1099 
986 
997 


1232 
1218 
1203 


1167 
1150 
1161 


1807 


6.11 x 


31 
31 


1997 


5.95 x 


34 


1029 


1934 


5.79 x 


10" 2 


- 


33 


31 


40 


40 


94 


- 


1012 


1008 


1218 


1203 


1892 


5.64 x 


10~ 2 


- 


33 


31 


40 


37 


92 


- 


1026 


1018 


1233 


1128 


1851 


5.48 x 


10~ 2 
10" 2 


- 


32 
32 


31 
31 


40 
39 


37 
37 


89 

87 


- 


1009 
1024 


1028 
1038 


1249 
1234 


1139 

1150 


1789 


5.33 x 


13 


737 


1749 


5.17 x 


10~ 2 


13 


32 


31 


39 


37 


86 


759 


1038 


1047 


1249 


1161 


1730 


5.02 x 


10~ 2 


14 


33 


32 


39 


37 


85 


843 


1085 


1090 


1267 


1171 


1711 


4.86 x 


10" 2 


15 


33 


32 


39 


37 


83 


929 


1100 


1099 


1283 


1182 


1672 


4.71 x 


10" 2 


16 


33 


32 


39 


37 


82 


1019 


1116 


1107 


1300 


1192 


1655 



TABLE IV: Proton and neutron number in the pasta phase region for several nuclear interactions. The horizontal lines represent 
the shape transitions, droplet-rod, rod-slab and slab-core. For NL3 only the droplet-core transition occurs, and for IU-FSU the 
droplet-rod transition occurs at the lowest density shown. 



1.2 p 
1.0- 

■■o.s- 



^0-6- 



0.4- 
0.2- 



NL3 

DD-ME2 

DD-MES 

NL3 o>p 

FSU 

IU-FSU 



0.0, 








FIG. 2: (Colors online) Surface tension as a function of the 
proton fraction x for all the models considered in this study 
except SLy4. 
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p 


- 3 ] 






Rw 


[fm] 










E/A [MeV] 






[fm 


NL3 DD-ME2 DD-ME5 NL3^ P 


FSU IU-FSU 


NL3 


DD-ME2 DD-ME5 NL3^ P 


FSU 


IU-FSU 


8.94 x 


10" 2 
10~ 2 








— 




- 








10.043 
10.000 




10.945 


8.88 x 


23.6 


10.922 


8.35 x 


10" 2 
10~ 2 








- 




16.9 

16.5 








9.614 
9.570 




10.693 


8.29 x 


11.7 


10.667 


7.82 x 


10 -2 


9.2 


12.0 


10.5 


11.5 


16.1 


15.4 


7.512 


9.351 


9.167 


9.217 


9.241 


10.448 


7.66 x 


10" 2 
10~ 2 


8.0 

8.8 


13.0 
11.0 


7.0 


11.7 


2.0 
10.2 


15.4 
15.4 


7.324 
7.141 


9.230 
9.110 


9.060 


9.100 
8.983 


9.130 
9.021 


10.373 


7.51 x 


10.7 


12.0 


8.951 


10.296 


7.35 x 


10~ 2 
10~ 2 
10~ 2 


9.3 

7.9 
8.9 


10.0 


10.8 
11.0 


14.1 
14.0 


11.8 
11.5 
11.5 


15.3 
15.5 
15.6 


6.961 
6.786 
6.614 


8.989 


8.842 
8.734 


8.865 
8.746 
8.627 


8.913 
8.802 
8.691 


10.217 


7.20 x 


11.7 

11.5 


8.868 
8.745 


10.137 


7.04 x 


13.1 


14.1 


8.626 


10.055 


6.88 x 


10 -2 
10~ 2 


8.4 
9.0 


11.7 

13.8 


13.1 
13.1 


14.1 
14.4 


11.6 
12.0 


15.8 
16.1 


6.445 
6.281 


8.622 


8.517 
8.409 


8.507 
8.386 


8.581 
8.470 


9.972 


6.73 x 


8.498 


9.886 


6.57 x 


10~ 2 


9.3 


13.7 


13.4 


14.5 


14.0 


16.3 


6.121 


8.373 


8.300 


8.264 


8.359 


9.798 


6.42 x 


10~ 2 


8.7 


13.7 


13.6 


14.7 


14.0 


16.6 


5.964 


8.248 


8.192 


8.142 


8.248 


9.708 


6.26 x 


10~ 2 
10~ 2 
10~ 2 


9.4 
9.3 
9.6 


14.0 
14.3 


13.9 


16.9 


14.3 
14.4 
14.6 


16.8 


5.812 
5.664 
5.519 


8.123 
7.998 


8.083 


8.019 
7.896 
7.771 


8.136 
8.024 
7.911 


9.616 


6.11 x 


15.8 
16.0 


17.0 
17.1 


19.1 
19.0 


7.974 
7.866 


9.522 


5.95 x 


16.2 


7.871 


9.423 


5.79 x 


10~ 2 


9.9 


16.3 


16.2 


17.3 


15.1 


19.1 


5.379 


7.745 


7.757 


7.646 


7.798 


9.322 


5.64 x 


10~ 2 


19.3 


16.5 


16.4 


17.5 


17.0 


19.1 


5.241 


7.617 


7.647 


7.520 


7.685 


9.219 


5.48 x 


10~ 2 


20.6 


16.5 


16.6 


17.8 


17.2 


19.1 


5.108 


7.490 


7.538 


7.394 


7.570 


9.112 


5.33 x 


10~ 2 


15.0 


16.8 


16.9 


17.9 


17.5 


19.1 


4.977 


7.362 


7.429 


7.266 


7.455 


9.003 


5.17 x 


10~ 2 


15.3 


17.0 


17.1 


18.1 


17.7 


19.3 


4.850 


7.234 


7.319 


7.138 


7.339 


8.890 


5.02 x 


10" 2 


16.0 


17.5 


17.5 


18.4 


17.9 


19.5 


4.726 


7.105 


7.208 


7.010 


7.223 


8.774 


4.86 x 


10" 2 


16.7 


17.7 


17.7 


18.7 


18.2 


19.6 


4.606 


6.976 


7.097 


6.880 


7.105 


8.654 


4.71 x 


10" 2 


17.4 


18.0 


17.9 


18.9 


18.4 


19.8 


4.491 


6.847 


6.986 


6.750 


6.987 


8.531 



TABLE V: Cell Wigner-Seitz radius and energy per baryon in the pasta phase region. The horizontal lines represent the shape 
transitions, droplet-rod, rod-slab and slab-core. For NL3 only the droplet-core transition occurs, and for IU-FSU the droplet-rod 
transition occurs at the lowest density shown. 
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p 


- 3 ] 








VN [M 


eV] 










X 








[fm 


NL3 


DD-ME2 DD-ME5 NL3„ P 


FSU 


IU-FSU 


NL3 


DD-ME2 DD-ME5 NL3^ P 


FSU 


IU-FSU 


8.94 x 


10" 


-2 








15.930 




13.711 








0.040 




0.045 


8.88 x 


10" 


-2 








15.837 




13.941 








0.040 




0.045 


8.35 x 


10" 


-2 








15.006 




13.752 








0.038 




0.046 


8.29 x 


10" 


-2 








15.049 




13.730 








0.038 




0.046 


7.82 x 


10" 


-2 


16.359 


14.759 


13.924 


14.419 


14.150 


13.529 


0.029 


0.037 


0.035 


0.037 


0.037 


0.047 


7.66 x 


10" 
10" 


-2 
-2 


15.800 
15.250 


14.518 
14.278 


13.672 
13.534 


14.207 
13.995 


13.876 
13.613 


13.460 
13.387 


0.028 
0.027 


0.037 
0.036 


0.034 


0.037 
0.036 


0.037 
0.036 


0.047 


7.51 x 


0.034 


0.047 


7.35 x 


10" 
10" 
10" 


-2 
-2 
-2 


14.712 
14.180 
13.661 


14.037 
13.867 
13.641 


13.310 
13.084 
12.891 


13.822 
13.613 
13.402 


13.473 
13.255 
13.034 


13.312 
13.235 
13.155 


0.027 
0.026 
0.025 


0.036 


0.033 
0.033 


0.036 
0.035 
0.035 


0.036 
0.036 
0.035 


0.047 


7.20 x 


0.036 
0.035 


0.047 


7.04 x 


0.033 


0.047 


6.88 x 


10" 
10" 


-2 

-2 


13.154 
12.654 


13.412 
13.216 


12.673 
12.454 


13.189 
12.976 


12.814 
12.595 


13.073 
12.989 


0.024 
0.023 


0.035 


0.032 
0.032 


0.035 
0.034 


0.035 
0.034 


0.047 


6.73 x 


0.034 


0.047 


6.57 x 


10" 


-2 


12.166 


12.987 


12.236 


12.759 


12.417 


12.902 


0.023 


0.034 


0.032 


0.034 


0.034 


0.047 


6.42 x 


10" 


-2 


11.687 


12.756 


12.019 


12.543 


12.204 


12.813 


0.022 


0.033 


0.031 


0.033 


0.034 


0.047 


6.26 x 


10" 
10" 
10" 


-2 
-2 
-2 


11.219 
10.763 
10.318 


12.521 
12.287 
12.064 


11.803 
11.597 
11.384 


12.338 
12.119 
11.896 


11.996 
11.785 
11.575 


12.718 
12.679 
12.584 


0.021 
0.020 
0.020 


0.033 
0.033 


0.031 


0.033 
0.033 
0.032 


0.033 
0.033 
0.033 


0.047 


6.11 x 


0.030 
0.030 


0.047 


5.95 x 


0.032 


0.047 


5.79 x 


10" 


-2 


9.884 


11.825 


11.172 


11.673 


11.368 


12.485 


0.019 


0.032 


0.030 


0.032 


0.032 


0.047 


5.64 x 


10" 


-2 


9.479 


11.585 


10.960 


11.447 


11.172 


12.382 


0.018 


0.031 


0.030 


0.031 


0.032 


0.047 


5.48 x 


10" 


-2 


9.087 


11.344 


10.750 


11.219 


10.966 


12.274 


0.018 


0.031 


0.029 


0.031 


0.031 


0.047 


5.33 x 


10" 


-2 


8.752 


11.102 


10.541 


10.990 


10.760 


12.162 


0.017 


0.030 


0.029 


0.031 


0.031 


0.047 


5.17 x 


10" 


-2 


8.370 


10.857 


10.333 


10.757 


10.555 


12.045 


0.017 


0.030 


0.029 


0.030 


0.031 


0.047 


5.02 x 


10" 


-2 


7.999 


10.612 


10.127 


10.524 


10.349 


11.922 


0.016 


0.030 


0.029 


0.030 


0.031 


0.047 


4.86 x 


10" 


-2 


7.637 


10.365 


9.922 


10.288 


10.144 


11.793 


0.016 


0.029 


0.028 


0.030 


0.030 


0.047 


4.71 x 


10" 


-2 


7.285 


10.118 


9.717 


10.050 


9.938 


11.659 


0.015 


0.029 


0.028 


0.029 


0.030 


0.047 



TABLE VI: Neutron chemical potential and proton fraction in the pasta phase region. The horizontal lines represent the shape 
transitions, droplet-rod, rod-slab and slab-core. For NL3 only the droplet-core transition occurs, and for IU-FSU the droplet-rod 
transition occurs at the lowest density shown. 



Transition 



droplet - rod 

rod - slab 

slab - homogeneous 



SLy4(HF) SLy4(TF) NL3 DD-ME2 DD-ME5 NL3^ P FSU IU-FSU 



0.042 
0.070 
0.080 



0.066 
0.078 
0.085 



0.0548* 



0.0611 


0.0626 


0.0642 


0.0580 


0.0471 


0.0688 


0.0720 


0.0751 


0.0673 


0.0626 


0.0735 


0.0766 


0.0835 


0.0751 


0.0894 



TABLE VII: Densities at the shape transition. The densities are give in fm ' 5 . 
found in [3]. For NL3, only the transition droplet-core is indicated. 



SLy4(HF) and SLy4(TF) refer to the results 
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FIG. 3: (Colors online) Neutron density at the cell center and 
border (a), cluster proton fraction at the cluster center (b) 
and neutron skin thickness (c) for all the models considered 
in this study. 
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FIG. 4: (Colors online) Proton number Z (a), neutron number N (b), Wigner-Seitz cell radius R ws (c), proton fraction x (d), 
the baryon number in the nucleus A (e) and the proton fraction in the nucleus Z/A (f) as a function of the density (p). HFB 
SLy4, HF N&V and CLD SLy4 refer respectively to the results found in [§] , [I[ and 0] . 
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FIG. 5: (Colors online) Proton number (a), neutron number 
(b) and Wigner-Seitz cell radius (c) in the pasta phase regions. 
CLD SLy4 refers to the results found in 



